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ABSTRACT 


As  annual  operating  efficiency  of  vented  heating  equipment  is  affected  by  burner 
fuel  and  combustion  air  modulation,  it  is  important  to  differentiate  between  the 
various  types  of  controls  in  determining  annual  energy  requirements.  Test  pro- 
cedures for  evaluating  annual  efficiency  have  already  been  developed  and  imple- 
mented by  the  Department  of  Energy  (DoE)  for  furnaces  with  single-stage  thermo- 
stat control.  A modified  test  procedure  is  necessary  to  account  for  operation 
with  fuel  modulation.  A revised  procedure  which  accommodates  two  types  of  fuel 
modulating  controls  has  recently  been  developed.  Tests  are  conducted  at  reduced 
and  maximum  firing  rates,  and  along  with  typical  derived  values,  from  a bin 
analysis  of  weather  data,  the  fraction  of  the  total  hours  for  each  operating 
mode  is  obtained  to  calculate  a weighted  annual  efficiency.  These  test  methods 
and  calculation  procedures  are  based  on  and  are  an  extension  to  the  current  DoE 
test  procedures  for  the  single-stage  type  of  thermostat  control  of  central  warm 
air  furnaces. 

By  using  the  procedures  developed  in  the  report,  the  energy  savings  impact  of 
fuel  modulating  controls  when  combined  with  the  use  of  modulated  combustion  air 
is  evaluated.  Energy  savings  from  6 percent  to  20  percent  were  determined  from 
the  increase  in  efficiency  with  both  fuel  and  combustion  air  modulation. 

Improved  efficiency  was  dependent  on  the  type  of  thermostat  control  and  the 
minimum-to-maximum  fueled  input;  i.e.,  turndown  ratio. 

Key  Words;  annual  efficiency;  household  heaters  and  furnace  test  procedures; 

hydraulic  thermostat  control;  modulating  control  gas-fueled; 
two-stage  thermostat. 
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INTRODUCTION 


As  annual  operating  efficiency  of  vented  heating  equipment  is  affected  by 
burner  operating  modes,  it  is  important  to  differentiate  between  the  various 
types  of  controls  in  determining  annual  energy  requirements.  Test  procedures 
for  evaluating  annual  efficiency  have  already  been  developed  and  implemented 
by  the  Department  of  Energy  (DoE)  for  heaters  with  single-stage  control  [3]*. 

I For  two-stage  and  step-modulating  units,  however,  a modified  test  procedure  is 
I necessary  to  account  for  operation  at  multiple  firing  rates,  A revised  proce- 
dure which  accommodates  modulating-type  units  has  recently  been  developed  and 
1 recommended  to  DoE  [1,  2].  A detailed  review  of  the  background  for  the  recom- 
mendations of  tests  for  modulating  heaters  in  that  procedure  is  the  subject  of 
this  report. 

jj  Modulating  controls  investigated  in  this  study  are  generally  used  with  vented 
household  heaters  other  than  furnaces.  However,  the  methods  of  test  and  the 
Ll  bin  analysis  procedures  developed  in  this  report  for  vented  heaters  are  also 
applicable  to  central  heating  systems  using  the  types  of  modulating  controls 
described  here.  The  only  adjustment  that  would  be  needed  to  these  procedures 
would  be  to  use  bin  data  based  on  the  population  of  central  furnaces  in 

j;  order  to  arrive  at  a national  average  annual  efficiency  specific  to  furnaces. 

jl 

1 

I Three  types  of  thermostatic  controls  are  commonly  used  on  automatically- 
i controlled  heating  equipment.  These  are; 

;! 

1 a.  single-stage  control,  which  cycles  the  heater  between  full  input  (100 
! percent  of  rating)  and  off; 

b.  two-stage  control,  which  operates  in  an  on-off  cycling  mode  at  either 

, maximum  input  or  at  some  reduced  input  which  can  be  as  low  as  20  percent 

of  maximum;  and 

c.  step-modulating  control,  which  steps  on  to  a low  input  and  then  either 

, cycles  off  and  on  at  the  low  input  if  the  heating  load  is  light,  or  grad- 

, ually  increases  the  heat  input  to  meet  any  higher  heating  load  that  can- 

' not  be  met  with  the  low  firing  rate. 

( 

I Figure  1 illustrates  the  burner  heat  output  rate  and  burner  on  periods  (shaded) 
area)  in  response  to  increasing  heating  loads  for  each  of  these  three  types  of 
controls. 

The  assumptions  made  in  previous  work  [4]  and  as  depicted  by  figure  1 that 
outdoor  temperatures  below  65°F  would  be  directly  proportional  to  heating  load 
is  valid  only  in  a broad  statistical  sense.  It  should  be  understood  that  proce- 
dures developed  here  for  calculating  annual  efficiency  as  well  as  current  DoE 
test  procedures  for  determining  annual  efficiency  of  the  conventional  type 
on-off  thermostatically  controlled  heaters  is  intended  for  comparison  purposes 


' Numbers  in  brackets  indicate  references  in  section  6. 

|i 

I 


I 

i 

I 
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between  equipment  tested  only  under  laboratory  test  conditions.  A comprehensive 
treatment  to  include  internal  heat  gains,  solar  radiation,  wind,  etc.,  would  be 
needed  in  order  to  be  valid  for  small  time  thermostatic  cycling  of  the  installed 
heating  system.  These  procedures  should  not  be  considered  as  being  predictive 
of  the  installed  annual  energy  use.  It  should  be  understood  that  a comprehen- 
sive treatment  of  the  building  and  furnace  interaction  to  predict  operating 
costs  after  installation  is  beyond  the  scope  of  this  work. 
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2 . STEP-MODULATING  CONTROLS 

2.1  MODE  OF  OPERATION 

One  representation  of  the  operating  mode  for  the  step-modulating  type  control 
is  its  response  to  room  temperature  (see  figure  2 where  room  temperature  and 
outlet  gas  pressure  supplied  to  the  burner  from  the  thermostat  control  are 
shown) . 

Under  normal  operating  conditions,  as  the  ambient  air  temperature  surrounding 
the  liquid-filled  thermostat  sensing  bulb  drops  to  point  A on  figure  2,  the 
control  opens  to  allow  minimum  fuel  flow  (point  B) . Uhen  room  air  temperature 
I rises  to  point  C,  the  valve  closes  (point  D) . Operations  on  the  ABCD  rectangle 
is  termed  the  "cycling"  mode, 

; If  more  severe  weather  is  encountered  while  the  valve  is  open  to  minimum  flow 
I (point  B),  room  temperature  may  continue  to  drop  until  point  E is  reached,  at 
which  time  any  further  drop  in  room  temperature  would  result  in  an  increasing 
; fuel  flow.  If  the  heating  demand  continues  to  increase  with  subsequent  lower 

i room  temperature,  fuel  flow  would  gradually  increase  to  maximum  flow  (from  E to 

j G) . In  another  situation,  if  the  heating  demand  should  just  equal  heater  out- 
; put  for  a prolonged  time,  the  manifold  could  remain  partially  open  (point  F or 
'1  any  other  point  between  low  fire  and  high  fire)  until  the  heating  demand  changes, 

i;  Operation  at  the  various  firing  rates  from  point  B to  G is  termed  the  "noncycling" 

' mode  or  modulating  mode. 

,1) 

2.2  ENERGY  CONSIDERATIONS 

In  developing  test  procedures  for  modulating  type  heaters,  operation  in  both 
, the  cycling  and  non-cycling  modes  must  be  addressed.  The  importance  of  addres- 
sing both  modes  separately  is  apparent  when  one  considers  the  difference  in 
efficiency  between  the  two  modes  as  illustrated  in  figure  3.  The  part-load 
^ efficiency  curve  of  figure  3 was  plotted  using  data  obtained  with  a room  heater 

j having  a maximum  rated  input  of  35,000  Btu/h  and  minimum  adjusted  input  of 

j 21,000  Btu/h  and  with  a balance  point  temperature  (T^)  of  36°F,  (See 
section  2.2.1  for  discussion  of  Tp.) 

I 

J Annual  efficiency  for  step-modulating  units  depends  upon  several  factors 
including: 

• percentage  of  the  heating  season  operating  time  in  each  of  the  two  modes, 

r 

• average  outdoor  temperature  for  each  mode, 

® average  steady-state  efficiency  over  the  various  firing  rates  in  the 

j;  non-cycling  mode, 

o infiltration  loss  in  the  non-cycling  mode,  and 

i' 

i|  • part-load  efficiency  in  the  cycling  mode. 


I 


The  first  two  factors  can  be  developed  by  considering  heater  sizing  and  through 
analysis  of  representative  weather  data.  The  remaining  factors  are  determined 
via  additional  test  measurements  and  calculations.  Each  of  the  factors  is 
described  in  detail  in  the  sections  that  follow. 

2.2.1  Percent  of  Operating  Time  in  Each  Mode 

Two  assumptions  provide  a basis  for  defining  the  percent  of  operating  time 
spent  in  each  of  the  two  modes:  (1)  at  an  outdoor  temperature  of  65°F,  heating 

requirements  no  longer  exist,  and  (2)  at  the  outdoor  design  temperature  for 
vented  room  heaters  of  15°F  (see  appendix  A for  development  of  this  value)  the 
heater  output  rate  needed  to  meet  the  heating  requirement  is  equal  to  the  maxi- 
mum heater  output  rate  divided  by  (1+a)  where  a is  the  heater  oversize  fraction. 
Here  output  rate  is  defined  as  the  product  of  input  rate  and  steady-state  effi- 
ciency. A relationship  between  heater  output  and  outdoor  temperature  can  then 
be  expressed: 


whe  re 


^out ,max/ ( 1+g)  _ 65-15 
Qout,min 


(1) 


Qout  max  rated  maximum  output  rate, 

Qout ,min  rated  minimum  output  rate. 

T^  is  the  outdoor  temperature  where  the  modulating  mode  begins  (see  figure  4). 
This  implies  that  heating  load  is  a linear  function  of  outdoor  temperature.  In 
these  test  procedures  a is  assigned  a value  of  zero.  Under  this  assignment,  at 
outdoor  design  temperature  of  15°F  the  heater  operates  at  maximum  firing  rate. 


A plot  of  Tj,  versus  the  ratio  of  minimum  to  maximum  output  is  shown  in  figure  5 
for  fixed  values  of  a.  The  cycling  mode  (as  shown  in  figure  4)  is  between  out- 
door temperature  T^  and  65° F.  The  modulating  mode  is  between  T(,  and  the  minimum 
outdoor  temperature,  which  for  household  heaters  is  typically  15°F,  Through 
analysis  of  typical  weather  data,  the  time  in  each  mode  can  then  be  defined  as 
a fraction  of  total  heating  season  degree  hours.  The  analysis  involves  a bin 
method  calculation  of  hours  in  various  temperature  ranges.  The  steps  in  the 
calculation  are  presented  in  the  column  headings  of  table  1.  See  column  (14) 
and  columns  (1)  through  (5)  of  table  1,  A graphical  example  is  also  shown  in 
figure  6. 

2.2.2  Average  Outdoor  Temperature  in  Each  Mode 


The  average  outdoor  temperature  in  the  cycling  and  non-cycliiig  modes  Tq^  and 
Toa'  > ^respectively  (as  shown  in  figures  7 and  8)  are  developed  from  representa- 
tive national  average  weather  data.  Tabulated  data  and  temperature  calculations] 
for  one  city  having  a 4200  degree  day  heating  season  are  shown  in  table  1.  The 
average  of  calculated  values  for  eight  cities  and  the  basis  for  figures  7 and  8 
are  presented  in  table  two.  The  eight  cities  are  listed  in  appendix,  table  A2 . 
These  eight  cities  all  fall  within  a range  of  +5  percent  of  the  national  average 
(heater  population  weighted)  of  4400  degree  days. 
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Referring  to  table  1,  calculation  of  Tq^  and  Tq^’  is  as  follows:  for  each 

temperature  range,  the  mean  temperature  (col.  2)  and  number  of  hours  spent  in 
the  range  (col.  3)  are  determined  from  a compilation  of  historial  weather  data 
[5].  From  this  the  number  of  degree  hours  in  each  range  (col.  5)  is  obtained. 
Toa  represents  the  average  outdoor  temperature  between  65°F  and  T^.  To  obtain 
Tq3,  col.  5 and  col.  3 are  therefore  summed  cumulatively  (upward  in  table  1) 
from  62°F  to  yielding  col.  6 and  col.  7,  respectively.  For  each  given  in 
col.  2,  col.  6 divided  by  col.  7 represents  the  corresponding  number  of  degrees 
below  65°F.  This  number  is  given  in  col.  8.  Toa  9)  is  65°F  minus 

col.  8.  The  same  approach  is  taken  in  calculating  Toa'  except  that  the  cumula- 
tive summation  is  from  15°F  to  T^.  rather  than  from  65°F  to  T^.  The  result  of 
these  calculations  is  a value  of  Toa  ^oa*  each  T^.  A plot  of  these 

parameters  for  the  values  in  table  2 is  shown  in  figure  7 . For  any  given 
heater,  equation  (1)  can  be  used  to  determine  Tq  , for  which  values  of  Toa 
Toa'  obtained  from  figure  7.  Alternatively,  equation  (1)  and  the  results 

of  the  previous  temperature  calculations  can  be  combined  to  eliminate  the  inter- 
mediate calculations  of  T^.  Toa  ^oa'  obtained  from  figure  8 where 

percent  of  time  in  cycling  modes,  Toa  ^oa'  (columns  5,  6a,  and  6b  of  table  2 
respectively)  are  plotted  vs  the  min/max  heater  output  ratio  (col.  7 of  table  2) 

2.2.3  Average  Steady-State  Efficiency 

Due  to  an  Increase  in  excess  air  at  below  maximum  firing  rates,  there  is  usually 
a drop  in  steady-state  efficiency  in  the  non-cycling  mode  as  fuel  input  rate 
is  reduced.  An  average  steady-state  efficiency  is,  therefore,  needed  to  repre- 
sent operation  in  this  mode.  The  average  steady-state  efficiency  is  calculated 
by  linear  interpolation  between  the  steady-state  efficiency  at  the  minimum 
rated  input,  which  corresponds  to  temperature,  T^,  and  the  steady-state  effi- 
ciency at  the  maximum  rated  input  for  the  heater,  which  corresponds  to  a 
temperature  of  15°F.  These  efficiencies  are  determined  in  accordance  with  the 
DoE  test  procedure  using  a=0.*  The  average  efficiency  (riss^  will  correspond  to 
the  efficiency  when  firing  at  temperture  Tq^'  and  is  given  by 


'ss 


'ss  ,max 


hss  ,minj 


Tc  - Toa’ 
T^  - 15 


+ n 


ss ,min 


(2) 


In  general,  the  arithmetic  mean  of  the  steady-state  efficiencies  may  be  used 
in  place  of  equation  (2)  with  a resulting  error  of  about  one  percentage  point. 


The  DoE  test  procedure  [3]  lists  a value  of  a=0.7  which  was  applicable  to 
central  furnaces,  and  was  originally  adopted  for  use  with  vented  heaters. 
Since  the  typical  oversize  factor  is  unknown  for  vented  heaters  the  use  of 
a=0.7  would  be  arbitrary.  In  order  to  reduce  the  complexity  of  the  calcu- 
lation procedures  for  modulating  type  heaters,  a value  of  a=o  is  applied 
for  all  calculated  efficiencies  reported  here. 
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2.2.4  Infiltration  Loss  in  the  Modulating  Mode 

On-period  infiltration  loss  (Lj  energy  efficiency  debit  for  heaters 

using  indoor  air  for  combustion  and  draft  hood  dilution.  The  loss  of  effi- 
ciency due  to  infiltration  of  outdoor  air  at  Toa'  (which  is  subsequently  heated 
to  an  indoor  temperature  of  typically  70°F)  is  subtracted  from  the  weighted 
average  steady-state  efficiency  calculated  by  the  above  procedure.  Average  on- 
period  infiltration  loss  is  calculated  in  the  same  manner  as  steady-state  effi- 
ciency. On-period  infiltration  losses  are  obtained  via  the  current  DoE  test 
procedure  [3]  and  the  average  infiltration  is  taken  to  be  the  arithmetic  mean 
of  the  values  determined  at  the  maximum  and  minimum  input  rates. 


3. 


TEST  PROCEDURE  AND  CALCULATION 


3.1  CURRENT  TEST  PROCEDURE 

The  current  test  procedure  for  determining  the  annual  efficiency  of  vented 
heaters  [ 3]  applies  only  to  units  equipped  with  single-stage  controls.  While 
not  directly  applicable  to  modulating-type  units,  only  minor  modifications  to 
the  procedure  are  necessary  to  accommodate  these  types  of  controls.  A cursory 
review  of  the  procedure  will,  therefore,  be  provided  here.  A more  detailed 
description  of  the  procedure  is  provided  in  references  [3,  4]. 

The  existing  procedure  involves  conducting  a steady-state  performance  test  on 
the  heater,  plus  measuring  flue  gas  temperatures  at  specific  times  during  the 
heat-up  and  cool-down  from  steady-state  conditions.  Specific  measurements  at 
steady-state  include  flue  and  stack  gas  temperature  and  CO2  concentration.  Two 
temperature  measurements  during  heat-up  and  three  during  cool-down  provide  the 
basis  for  exponential  approximations  for  the  temperature-time  profiles. 

In  addition  to  the  above  experimental  data,  a knowledge  of  several  factors 
describing  mass  flow  rates  through  the  flue  and  stack  during  on-  and  off-periods 
is  required.  Values  for  these  factors  may  be  either  measured  or  assigned 
according  to  the  type  of  unit  under  test. 

Based  upon  the  experimental  data  and  assigned  system  factors,  thermal  losses 
associated  with  heater  operation  are  determined.  These  losses  expressed  as  a 
percentage  of  the  fuel  input  rate  are: 

Ll  = Latent  Heat  Loss  due  to  the  presence  of  uncondensed  water  vapor  in 
the  flue  gas. 

^S,ON  “ On-cycle  Sensible  Heat  Loss  due  to  the  venting  of  combustion  products 
and  excess  air  at  a temperature  above  room  temperature. 

Ll,ON  ~ On-cycle  Infiltration  Loss  due  to  heating  the  on-period  combustion 
and  draft  control  air  from  outdoor  temperature  (Tqa)  to  room 
temperature. 

Ls,off  “ Off-period  Sensible  Heat  Loss  due  to  heating  the  off-period  room  air 
discharged  at  a temperature  in  excess  of  the  room  temperature. 

^I,OFF  “ Off-period  Infiltration  Heat  Loss  due  to  heating  the  off-period  room 
air  loss  from  outdoor  temperature  to  room  temperature. 

Steady-state  efficiency  is  given  in  terms  of  latent  heat  loss,  Lj^,  and  on-perlod 
sensible  heat  loss  at  steady-state,  Lg  qn  ss*  Part-load  efficiency,  is 

expressed  as  a function  of  the  five  losses: 

riu  = 100  - Ll  - — (Ls,0N  + 1^S,0FF  + Ll,0N  + Li,0FF)  (3) 

^on  X ’^off 
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typical  on-period  time  - minutes, 
typical  off-period  time  - minutes. 

Pilot  fraction  (as  fraction  of  total  fuel  input  rate). 


^on 


^of  f 


PF 


The  part-load  efficiency  is  then  combined  with  the  steady-state  efficiency 
(pgs)  and  pilot  fraction  to  yield  the  annual  fuel  utilization  efficiency. 


EFFYA  = 


^ss  ^ ^u  ^ 

Pgg  X 4400  + 2.5  X riu  X PF  X 4600 


(4) 


where  the  additional  parameters  are: 

4400  = average  annual  degree  days  for  vented  heaters  (see  appendix  A) 

4600  = average  non-heating  season  hours  per  year  that  all  the  energy 

to  the  pilot  is  assumed  wasted. 


3 . 2 MEASUREMENTS  AND  CALCULATIONS  FOR  HEATERS  WITH  STEP-MODULATING  CONTROLS 


For  single-stage  heaters  there  is  only  one  fuel  input  rate  at  which  measure- 
ments can  be  made  (the  maximum  input  rate).  Modulating- type  heaters,  however, 
operate  both  at  reduced  firing  rate  during  the  cycling  mode  and  anywhere  between 
the  reduced  and  maximum  rates  during  the  non— cycling  mode.  Accordingly,  test 
measurements  must  be  made  at  the  maximum  input  rate  as  well  as  the  reduced  set- 
ting. Since  the  heater  with  step-modulating  control  cycles  on-off  only  at  the 
reduced  input  rate,  all  part-load  cycling  losses  must  be  determined  at  the 
reduced  input.  In  addition,  to  account  for  operation  in  the  non-cycling  mode, 
as  described  in  sections  2.2.3  and  2.2.4,  Lg^oN  ^I,0N  be  determined  at 

maximum  fire  as  well.  It  is  not  necessary  to  know  either  of  the  off-cycle 
losses  at  high  fire  (Lg  qfF  OFF^  because  no  off  period  loss  occurs  at 

that  input  rate  with  the  step-moduiating-type  control. 

The  recommended  procedure  for  calculating  the  annual  fuel  utilization 
efficiency  for  modulating  type  heaters  essentially  involves  determination  of 
losses  and  efficiencies  for  each  of  the  operating  modes  — cycling  and  non- 
cycling and  weighting  of  these  parameters  according  to  the  fraction  of  time 
spent  in  each  mode. 

The  procedure  for  calculating  annual  efficiency  for  the  step-modulating  heaters 
is  summarized  below  in  nine  steps.  The  detailed  procedure  consisting  of  20 
steps  is  included  in  appendix  B,  along  with  a sample  calculation.  A computer 
program  (FBVH)  for  conducting  the  complete  calculation  applicable  to  furnaces, 
boilers,  and  vented  heaters  is  included  in  appendix  C. 

(a)  Determine  minimum  and  maximum  heater  outputs  from  the  minimum  and  maximum 
input  rates  and  measured  steady-state  efficiency  at  these  two  input 
rates.  (For  details  see  steps  1-7  of  appendix  B.) 
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(b)  Using  Tnlniraum/ maximum  heater  output  from  step  (a) , determine  the  percent 
of  heating  season  in  each  mode.  Find  the  percent  of  time  in  cycling  mode 
from  figure  8.  Percent  of  non-cycling  mode  is  100  percent  less  the 
percentage  of  cycling  mode. 

(c)  Determine  the  average  outdoor  temperature  in  the  cycling  mode,  Tq^,  and  in 
the  non-cycling  mode,  Tq^' , from  figure  8 at  the  point  corresponding  to  the 
min/max  output  determined  from  (a). 

(d)  Determine  part-load  efficiency  in  the  cycling  mode  using  the  value  of  Toa 
and  the  prescribed  DoE  test  procedure  [3]. 

(e)  Determine  the  average  steady-state  efficiency  in  the  non-cycling  mode. 

Use  the  average  of  steady-state  efficiency  measured  at  the  maximum  and 
minimum  input  rates  from  step  (a). 

(f)  Determine  part-load  efficiency  in  the  non-cycling  mode  by  subtracting 

infiltration  loss,  qN»  the  average  steady-state  efficiency  in 

step  (e).  Use  Tq^* , to  calculate  infiltration  loss  in  the  non-cycling  mode. 

(g)  Determine  weighted  average  steady-state  efficiency  using  steady-state 
efficiencies  at  the  minimum  input  rate  from  step  (a)  and  the  average  deter- 
mined for  the  non-cycling  mode  from  step  (e).  Weight  each  by  the  percent 
of  heating  season  in  each  mode  (determined  from  step  (b)). 

(h)  Determine  average  part-load  efficiency  for  the  heating  season  using  the 
corresponding  part-load  efficiencies  for  the  cycling  and  non-cycling  modes 
and  the  percent  of  heating  season  in  each  mode. 

(i)  Determine  annual  efficiency  using  average  part-load  efficiency  and  average 
steady-state  efficiency. 

3.3  CALCULATIONS  FOR  HEATERS  WITH  TWO-STAGE  CONTROLS 


This  type  of  control  cycles  the  burner  at  low  fire  for  outside  air  temperatures 
between  65°F  and  T^  and  at  higher  fire  for  temperatures  T^  and  below  (see  figure 
1).  Evaluation  of  the  annual  fuel  efficiency  of  heaters  with  two-stage  controls 
involves  the  same  test  measurements  required  for  step-modulating  type  units,  and 
only  slightly  different  calculations.  The  calculations  are  described  in  detail 
in  appendix  B. 
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IMPACT  OF  MODULATING  CONTROLS  ON  ENERGY  CONSUMPTION 


4.1  PRESENT  DESIGN 

Figure  3 shows  that  modulating  controls  will  indeed  affect  annual  efficiency 
as  evidenced  by  the  variation  in  part  load  efficiency  with  heating  load.  As 
shown,  efficiency  improves  as  the  heating  load  is  increased.  This  is  due  both 
to  a reduction  in  heater  cycling  (see  Modulating  Mode  of  figure  4)  and  there- 
fore, off-cycle  losses,  and  also  in  the  higher  steady-state  efficiency  at  the 
higher  full  input  rates.  The  maximum  efficiency  of  currently  designed  heaters 
is  realized  in  the  non-cycling  mode. 

4.2  HEATERS  EQUIPPED  WITH  THERMAL  DAMPERS 

Measured  on-period  and  off-period  stack  losses  of  a room  heater  with  and  without 
a thermal  stack  damper  installed  were  measured.  Detailed  test  results  and  data 
are  in  tables  3 & 4.  Data  for  damper  A is  presented  graphically  in  figure  9. 

See  [6]  for  additional  information  including  description  of  tracer  gas  test 
method  use  to  measure  the  off  period  losses.  The  shaded  area  represents  energy 
savings  due  to  the  damper.  Figure  9 shows  for  that  heater,  as  fuel  input  rate 
is  reduced  from  maximum  to  lower  input  rates,  energy  savings  due  to  the  damper 
continuously  increase.  The  significance  of  these  findings  is  that  a test  con- 
ducted per  the  DoE  test  procedure  [3]  only  at  the  maximum  input  rate  does  not 
reflect  the  potential  energy  savings  of  the  thermal  type  stack  damper  when  it 
is  installed  on  a heater  having  a step-modulating  or  two-stage  thermostat 
control . 

A summary  of  energy  savings  calculated  using  the  calculation  procedure  of 
appendix  C for  three  heaters  with  three  different  models  of  thermal  dampers  is 
shown  in  table  5.  Results  presented  in  table  5 show  that  energy  savings  is 
dependent  on  type  of  heater  and  stack  damper  as  well  as  type  of  thermostat  used, 
and  the  reduced  fuel  input  rate  of  the  step-modulating  type  thermostat. 

4.3  HEATERS  WITH  REDUCED  COMBUSTION  AIR  AT  REDUCED  HEAT  INPUT  RATES 


Although  heaters  tested  have  shown  a reduced  part-load  efficiency  at  reduced 
input  rates,  it  is  possible  that  heaters  can  be  designed  to  use  less  excess  air 
at  the  reduced  input  rate,  thereby,  actually  increasing  part-load  cycling  mode 
efficiency  at  the  lower  rate  to  above  the  non-cycling  mode  operating  efficiency. 
This  is  possible  because  at  reduced  input  rates  the  ratio  of  heat  transfer  sur- 
face area  to  combustion  products  mass  flow  rate  is  greater.  Test  procedures 
outlined  here  will  allow  credit  for  any  such  innovative  designs  that  result  in 
reduced  excess  combustion  air  at  reduced  fuel  input  rate. 

The  potential  for  energy  savings  with  reduced  excess  air  has  been  demonstrated 
in  the  laboratory  by  intentionally  reducing  the  excess  air  at  reduced  input 
rate.  The  excess  combustion  air  values  at  a minimum  input  are  presented  in 
table  6 for  a heater  in  "as  received"  condition  and  after  reducing  excess  com- 
bustion air.  Excess  air  was  reduced  in  these  test  by  using  a flue  baffle  placed 
in  the  exit  of  the  heat  exchanger.  Baffling  of  the  heat  exchanger  was  limited 
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in  order  not  to  increase  carbon  monoxide  in  the  flue  gases  above  the  amounts 
found  prior  to  any  restriction  of  the  flue.  Data  obtained  at  the  maximum  firing 
rate  and  at  the  minimum  input  rate  with  restricted  excess  combustion  air  is 
shown  in  table  6. 

The  effect  of  reduced  excess  air  on  efficiency  can  be  quite  significant.  As 
shown  in  table  6,  energy  savings  of  from  6 to  22  percent  were  calculated. 
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5. 


CONCLUSIONS 


A procedure  for  evaluating  the  annual  operating  efficiency  teen^devel- 

eaulpment  with  a step-modulating  thermostat  or  two-stage  control  has  been  devel 
oped.  This  procedure  is  essentially  the  same  as  the  existing 

dLe  for  venLd  heaters,  but  calls  for  tests  and  calculations  to  performed 
at  maximum  as  well  as  mlnumum  firing  rates.  Analyses  of  weather  data  typl 

cal  cities,  and  assumptions  concerning  heater  siring,  provide  a means  of  com 
blnlng  the  cycling  and  non-cycling  efficiencies  to  yield  a weighted  annual 

fuel  efficiency. 

This  rated  annual  efficiency  is  considered  to  be  reflective  of  differences 

Inherent  with  cycling  and  non-cycling  modes  of  operation  which  apply  with 
mLulatlng-type^ont?ols.  Since  procedures  developed  here 

current  DoE  test  procedures,  these  procedures  are  believed  to  be  “seful  o 
cLparing  the  annLl  efficiency  of  heaters  with  modulating  type 
conventional  on-of  thermostat  control  used  with  some  types  of 

The  bln  calculation  procedures  developed  here  are  also  expected  to  be  applicable 
w certain  furnaces  equipped  with  the  types  of  fuel  modulating  controls  Investi- 
gated in  this  study.  Predicted  annual  energy  with  any  of  these  ^ 

should  be  considered  valid  only  for  comparison  purposes  between  equipment  teste 
under  laboratory  conditions.  These  procedures  as  well  as  the  current  D 
procedures  should  not  be  considered  predictive  of  installed 
costs.  Comprehensive  treatment  of  the  building  and  furnace 

are  unique  to  each  installation  would  be  needed  for  that  prediction  and  Is 
beyond  the  scope  of  these  procedures. 

Comparison  of  heater  efficiencies  with  single-stage  and  with 

controls  Indicates  that  substantial  Improvements  can  be  made.  ® j 

combustion  air  to  fuel  ratio  with  reduced  Input  rate  is  a promlstng  means  of 
reducing  energy  consumption  for  heaters  with  step-modulating  and  two  stage 

controls. 
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Table  1.  Development  of  outdoor  temperatures  and  time  In  each  mode  from  weather  data  for  one  city  (Bishop,  California) 
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NOTE:  All  calculated  values  in  this  table  are  based  on  no  oversizing  (a  = 0) 

Each  line  of  Col,  6 represents  the  average  from  eight  cities  with  individually 
calculated  values  from  specific  data  for  each  city. 


Table  3.  Summary  of  stack  losses  and  efficiencies  for  a room  heater  rated  45,000  Btu/h  (13,185W) 
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Table  5.  Effect  of  type  of  thermostat  control  on  energy  saving  potential  of  thermal  dampers 
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NOTE:  For  additional  details  on  test  results  and  calculation  procedures  for 

heaters  equipped  with  thermal  dampers,  see  Ref  [6]. 


Table  6.  Effect  of  reduced  excess  combustion  air  on  efficiency  of  a 35,000  Btu/h  room  heater  (10,255W) 
equipped  with  (1)  step-modulating  thermostat,  or  (2)  two-stage  thermostat 
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HEATING  LOAD—  (outdoor  temperature) 

Figure  1.  Typical  operating  modes  for  heaters  equipped  with  single-stage, 
two-stage,  and  step-modulating  controls 


Figure  2 


ROOM  TEMPERATURE 
OR  SENSING  BULB  TEMP. 


Outlet  gas  pressure  from  a step-modulating  thermostat  control 
in  response  to  changing  room  temperature 


21 


o 


22 


lodulating  control 


Heater  Output  B 
Design  Heater  Output  A 


) 


O 

h" 

<1 


o 


Q> 


O 

>» 

o 

— 

o 

o 


■XBIAI  |0  o/^ 

e\B\i  indjno  lean 


o 


“O 

cc 

o 

O) 

£ 

■ MM 

CO 

0 


Q 

h" 

d 


O 

H 

<1 


o 

LO 

CD 


23 


Figure  4.  Burner  output  rate  versus  heating  load  and  associated  outdoor 
temperatures  for  a heater  with  step-modulating  control 


OUTDOOR  TEMPERATURE 
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Figure  5. 


Balance  point  temperature  as  a function  of  heater  output  ratio 
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Cycling  Mode  Fraction  = 2530  -f-  4400 

or  58% 

Non-Cycling  Mode  = 100-58  - 42% 
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Figure  6.  Example  of  bin  method  used  to  determine  percentage  of  time 
in  cycling  mode 
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AVERAGE  OUTDOOR  TEMPERATURE,  IN  STEADY  STATE  AND  CYCLING  MODES 


Figure  7.  Average  outdoor  temperature  in  cycling  and  non-cycling  modes 
versus  balance  point  temperature  (T^^) 
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Figure  8. 


Average  outdoor  temperature  and  percent  of  time  in  cycling 
mode  as  a function  of  heater  output  ratio 
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Figure  9.  The  effect  of  burner  fuel  input  rate  on  the  part-load  loss  of 
room  heater  with  and  without  a thermal  damper 
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Appendix  A 


Determination  of  Average  Heating  Degree  Days  and  Design 
Temperature  for  Vented  Room  Heaters 

The  97  1/2  winter  design  temperature  is  that  dry-bulb  temperature  which  will 
probably  be  exceeded  97  1/2  percent  of  the  time  for  winter  months  (December, 
January,  February).  Through  consideration  of  the  geographic  distribution  of  l| 

vented  room  heaters  throughout  the  U.S.,  combined  with  analyses  of  weather 
data,  a single  value  or  design  temperature,  representative  of  a "National  ;j 

Average"  condition,  can  be  obtained.  This  value  will  be  different  than  the  ! 

corresponding  value  of  5°F  for  furnaces  used  in  [4] , since  the  geographic  ; 

distribution  of  furnaces  and  room  heaters  is  not  the  same. 

Based  on  1970  census  data,  the  total  number  of  room  heaters  installed  in  each 
of  the  nine  geographic  regions  of  the  United  States  is  presented  in  table  A1 . j 

Also  presented  in  table  A1  are  representative  average  number  of  heating  season 
degree  days  for  each  of  the  nine  regions.  Based  on  these  two  pieces  of  informa- 
tion, a heater  population  weighted  average  of  space  heating  degree  days  for 
room  heaters  can  be  obtained.  Results  of  the  average  are  presented  at  the 
bottom  line  of  column  (4)  (4438  is  rounded  to  4400  degree-days,  as  the  weighted 
national  average) . Although  these  data  are  based  on  the  1970  census  now  11  years  ' i 
old,  results  are  expected  to  be  little  affected  if  more  recent  census  data  were 
to  be  used  since  most  vented  heaters  sold  are  replacement  units  for  existing 
installations . 

The  outdoor  design  temperature  can  be  approximated  as  the  average  of  the 
design  temperatures  of  a representative  number  of  cities  having  heating  seasons 
of  about  4200  to  4600  degree-days.  The  heating  degree-days  and  outdoor  design 
temperatures  for  13  cities  falling  in  this  range  are  presented  in  table  A2. 

The  average  design  temperature  for  the  group  is  15.7°F. 
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Table  A-1 


Population  Weighted  Average  Space  Heating  Degree  Days 


(1) 

(2) 

(3) 

(4) 

Census 

Region 

Population 
of  Vented 
Heaters 

% of 
Total 

Regional 
Average 
Degree  Days 

Population 
Weighted 
Degree  Days 

Pacific 

3,672,950 

26.6 

3820 

1016 

Mountain 

649,618 

4.7 

5780 

271 

East  So.  Cntrl . 

1,007,956 

7.3 

3795 

277 

So.  Atl. 

2,683,320 

19.5 

3380 

659 

W.  No.  Cntrl. 

1,179,377 

8.6 

6825 

587 

E.  No.  Cntrl. 

1,640,383 

11.9 

6610 

787 

Mid.  Atl. 

800,859 

5.8 

6100 

354 

New  Eng. 

369,041 

2.7 

6820 

184 

W.  So.  Cntrl. 

1,783,762 

12.9 

2350 

303 

Total 

13,787,266 

4438 
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Table  A-2  Design  Temperatures  from  97  1/2%  Column  of  ASHRAE 
Handbook  of  Fundamentals  (1972) 

DEGREE 


CITY 

DESIGN  TEMP 
°F 

DAYS 

Charlestown,  SC 

14 

4476 

Huntington,  WV 

14 

4446 

Philadelphia,  PA 

15 

4486 

St.  Louis,  MO 

8 

4484 

Albuquerque,  NM 

17 

4348 

Evansville,  IN 

10 

4435 

Bishop,  CA 

15 

4275 

Washington,  DC 

15 

4551 

Sandberg,  CA 

25 

4209 

Eureka,  CA 

35 

4643 

Prescott,  AZ 

15 

4654 

Dover,  DE 

12 

4660 

Louisville,  KY 

9 

4620 

Avg. 

15.7 

APPENDIX  B 


Detailed  Stepwise  Procedure  for  Calculation  of  Annual 
Efficiency  for  Heaters  Equipped  with  Step-Modulating 
or  Two-Stage  Thermostat  Controls 
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Detailed  Stepwise  Procedure  for  Calculation  of  Annual 
Efficiency  for  Heaters  Equipped  with  Step-Modulating 
or  Two-Stage  Thermostat  Controls 


The  first  11  steps  deal  with  the  calculation  of  minimum  and  maximum  heat 
outputs,  and  fraction  of  time  and  outdoor  temperature  in  each  mode.  Each 
of  the  following  steps  refers  to  a similarly  numbered  block  in  figure  B-1. 


1 . Enter  the  maximum  fuel  input  rate  (including  fuel  supply  to 
pilot  flame) . 


2. 

3. 


Enter  the  steady  state  efficiency  at  the  maximum  input  rate 
(determined  in  accordance  with  the  DoE  test  procedure).* 


Calculate  and  enter  the  maximum  heat  output  rate 


Step  1 X Step  2 
100 


4.  Enter  the  minimum  (reduced)  fuel  input  rate  (including  fuel  supply 
to  pilot  flame) . 


5. 

6. 


Enter  the  steady  state  efficiency  at  reduced  fuel  input  rate 
(determined  in  accordance  with  the  DoE  test  procedure).* 


Calculate  and  enter  the  minimum  heat  output  rate 


Step  4 X Step  5 
100 


7.  Calculate  and  enter  ratio  of  minimum  to  maximum  output 


RATIO 


Q 

out ,min 

0 _ 

out ,max 


or 


Step  6 
Step  3 


8.  Read  the  percentage  of  time  in  the  low  fire  cycling  mode  fi^om  the 

graph  of  figure  7 for  the  minimum/maximum  ratio  corresponding  to  the 
ratio  determined  in  step  (7) . 

9.  Calculate  percentage  of  time  in  the  non-cycling  mode  (Rj^^)  . 

R^.  . 100  - (Step  8) 

10.  Read  the  average  outdoor  temperatures  for  the  cycling  and  non-cycling 

11.  mode,  T and  Tq^'  respectively,  from  figure  7 at  the  point  corresponding 
to  min/max  ratio  determined  in  step  (7) . 


In  the  next  four  steps,  the  average  steady  state  efficiency  and  on-cycle 
infiltration  losses  for  the  non-cycling  mode  are  determined. 


*References  1 and  3. 
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12.  Calculate  and  enter  average  steady  state  efficiency  (%)  for  the  non- 
cycling mode, 


13.  Enter  the  on-perlod  infiltration  loss  (Lj  0N,niin)  at  the  minimum  input 

rate  (determined  in  accordance  with  the  £>oE  test  procedure  for  an 
outdoor  temperature  equal  to  (Step  11), 

14.  Enter  the  on-cycle  infiltration  loss  at  the  maximum  input  rate 

(determined  in  accordance  with  the  DOE  test  procedure  for  an  outdoor 

temperature  equal  to  T , (Step  11) . 

oa 

15.  Calculate  and  enter  the  average  on-cycle  infiltration  loss  (%)  for  the 
non-cycling  mode 

, _ ON, min  ON, max  Step  13  + Step  14 

h,ON,avg  2 2 

The  part  load  efficiency,  n,j,  for  both  the  cycling  and  non-cycling  modes  are 
determined  next,  and  the  weighted  average  efficiencies  are  then  calculated. 


16.  Calculate  and  enter  the  part  load  efficiency  for  the  non-cycling 
mode 


n 


u,hi 


”ss,hi  (Step  12)  - I-I, 


ON, 


(Step  15) 
av 


17.  Calculate  and  enter  the  percentage  part  load  efficiency  (%)  for  the 
iwii^-cycling  mode : 

20 

^u,lo  " ■ t^^  + PFxt^pp  ^^S,0N  h,0N  S,0FF  h,0FF^ 

where  the  thermal  losses  are  determined  in  accordance  with  the  DOE 
nrocedure  for  the  minimum  firing  rate,  at  outdoor  temperature  equal 
to  T (Col  10),  and  t^j^  = tQpj.  = 20  min.  Alternatively ,Lg  Qpp  and 
Li  qpj,  can  be  measured  using  the  tracer  gas  technique  [6].  PF  is  the  pilot 
fraction  (the  ratio  of  pilot  input  Qp  to  total  fuel  input  . 

18.  Calculate  and  enter  the  weighted  average  part  load  efficiency  for  the 
cycling  and  non-cycling  modes 


u,wt 
or  Step 


R,  X n , + R,  , X n 

lo  u,lo  ni 

8 X Step  17  + Step  9, 


u ,hi 
Step 


16 


Appendix  B (cont) 
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19.  Calculate  and  enter  the  weighted  average  steady  state  efficiency  for 
the  cycling  and  non-cycling  inodes. 

"ss.wt  ' "lo  ==  "ss.mln  *hl  "ss,hl  ^ ^tep  5 + Step  9 X Step  12 

The  annual  fuel  efficiency  for  modulating  type  heaters  is  finally  calculated 
using  the  weighted  efficiencies  from  steps  18  and  19. 


20.  EFFYA  = 


ssjWt  X n 


U,wt 


X 4400 


n X 4400  + 2.5  n X PF  X 4600 

SSjWt  u,wt 


Procedure  for  Calculating  Annual  Efficiency  for  Heaters  with  Two-Stage  Controls 


The  calculations  listed  above  for  heaters  with  step-modulating  controls  will, 
with  slight  modification,  accommodate  heaters  with  two-state  controls.  The 
modifications  are  as  follows: 

o replace  step  (12)  with  the  equality:  n u*  ~ ^ 

s s y m ss^ni3x 

o delete  steps  (13)  through  (15) , and 

o replace  step  (3fT}  with  a calculation  similar  to  that  of  step  (JPfr) 

but  for  the  maximum  input  rate  and  an  outdoor  temperature  equal  to 


The  annual  fuel  efficiency  given  by  step  (20)  will  then  be  for  a heater  with 
two-stage  contol. 


Appendix  B (cont) 
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For  Step  Modulating  Type  Heaters  (Operating  in  a Cycling  Mode  at 
Reduced  Input  Rate  (Low  Fire)  and  Non-Cycling  Mode  at 
Intermediate  Input  Rates  Between  Low  and  Maximum  Rated  Input 
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Efficiency  and  Energy  Consumption  of  Vented  Heating  Equipment 
With  Single-Stage,  Two-Stage,  and  Step-Modulating  Thermostats 
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Computer  Program  for  the  Evaluation  of  Annual  Fuel  Efficiency  of  Heating 
Equipment  With  Single-Stage,  Two-Stage,  and  Step-Modulating  Control 


Program  Description 

The  computer  program  to  be  described  here  is  based  on  the  original  program 
NBSFBS5  developed  in  1978  by  Kelly,  et  al  [2].  While  maintaining  the  previous 
structure  of  NBSFBS5,  the  present  program  handles  heating  equipment  with  single- 
stage,  two-stage,  and  step-modulating  input,  and  permits  the  use  of  tracer  gas 
measurements  in  determining  off-cycle  losses  (see  Reference  ) . The  computer 
program  consists  of  a main  program,  entitled  FBVH,  and  four  subroutines:  SENLOS,* 
WEIGHT,  OFFLOS,  and  FUNT4.  It  is  written  in  ASCII  Fortran  and  should  be 
compatible  with  most  Fortran  V processors  with  a minimum  of  change. 

The  main  program  FBVH,  as  its  predecessor,  chooses  calculation  paths,  makes 
appropriate  calls  to  subroutines,  and  performs  calculation  steps  based  on  input 
information.  The  theoretical  calculations  and  the  order  in  which  the  calculations 
proceeds  are  essentially  the  same  as  in  NBSFBS5.  For  single-stage  controls,  a 
straight-through  path,  identical  to  that  in  NBSFBS5,  is  taken.  For  step-modulating 
and  two-stage  controls,  outdoor  temperatures  and  time  in  each  mode  are  determined 
via  a call  to  t,jeIGHT,  and  two  passes  through  the  calculation  stream  are  then 
taken.  In  the  first  pass,  input  data  for  the  maximum  firing  rate  and  outdoor 
temperature  T^^,  are  used.  In  the  second  pass,  required  parameters  are  saved, 
and  data  for  the  minimum  firing  rate  and  outdoor  temperature  To^f  are  used. 

Weighted  efficiencies  are  calculated  in  accordance  with  section  3.2. 

Subroutine  SENLOS  determines  values  for  steady  state:  latent  heat  loss 

sensible  heat  loss  Lg  and  steady  state  efficiency  Tgg  as  well  as  other 

parameters.  This  SENiOS ’version  contains  the  stack  to  flue  mass  flow  ratio  S/F 
(SFR)  computed  from  stoichiometric  relationships  and  is  used  in  place  of  the 
input  value  of  SFR  for  all  calculations  (the  computed  value  was  previously  used 
only  if  it  was  greater  than  the  input  value  of  SFR) . 


Subroutine  WEIGHT  assigns  values  for  average  outdoor  temperatures  in  the  cycling 
and  non  cycling  modes  (Toa  and  Toa'  respectively) , and  the  fraction  of  time  in  the 
cycling  mode.  WEIGHT  is  called  only  for  step-modulating  and  two-stage  heaters. 

The  subroutine  OFFLOS  calculates  off-cycle  sensible  heat  and  infiltration  losses 
based  on  the  tracer  gas  technique  described  in  Reference  6,  For  this  calculation, 
which  is  optional,  twenty  measurements  of  tracer  concentration  in  stack  gas  are 
required  at  one-minute  intervals,  along  with  the  initial  tracer  concentration,  and 
tracer  gas  flow  rate  and  temperature. 

Subroutine  FUNT  A calculates  functions  required  for  the  analytical  determination 
of  off-cycle  stack  losses  (the  current  DOE  and  default  means  of  calculating  off- 
cycle  losses) . FUNT  4 is  unchanged  from  the  NBSFBS5  version. 


*SENLOS  used  in  this  procedure  is  modified  from  another  form  of  the 
SENLOS  given  in  Reference  2 with  NBSFBS5. 
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A flow  chart  of  FBVH  including  corresponding  statement  numbers  from  the 
program  listing,  is  shown  in  figure  Cl.  This  figure  does  not  include  all 
program  details  but  only  those  essential  differences  between  FBVH  and  NBSFBS5. 
The  majority  of  the  theoretical  calculations  are  performed  in  statements  132 
through  256  and  are  described  in  detail  in  Reference  k . The  complete  program 
listing,  with  subroutines,  is  shown  in  figure  C2. 

Input  Data  Format 

The  input  data  file  or  element  for  FBVH  must  be  consistent  with  the  format 
shown  in  figure  C3.  A sample  data  form  is  presented  as  figure  C4  for  comparison 
It  should  be  noted  that  when  ICTRL  and  LOSOPT  are  selected  such  that  some  lines 
in  figure  C3  are  omitted  in  the  data  file,  blank  lines  should  not  be  inserted; 
the  line  numbers  in  the  data  file  will  therefore  differ  from  those  in  figure  C3 
in  certain  cases. 

Sample  Input  and  Output 

Sample  sets  of  computer  program  output  for  a heater  with  a step-modulating  and 
with  a two-stage  control  are  presented  in  figures  C4  and  C5  respectively. 


Appendix  C (cont) 
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Figure  Cl.  Flow  chart  of  main  program 
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Figure  Cl  (cont) 
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Figure  C2.  Main  program  and  subroutine  listings;  main  program  FBVH 


HTRS*B0IL6A( 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 
17 


1 ) .FBVH(3) 

C PROGRAM  FOR  THE  EVALUATION  OF  FURNACE/BOILER/VENTED  HEATER 

C SYSTEMS  WITH  SINGLE.  TWO-STAGE.  OR  STEP -MODU LATED  INPUT. 

C THIS  PROGRAM  IS  BASED  ON  THE  ORIGINAL  PROGRAM  NBSFBVH 

C (BOILS)  DEVELOPED  IN  1978  BY  KELLY . CHI . AND  KUKLEWICZ. 

C WHILE  MAINTAINING  THE  STRUCTURE  OF  B0IL6.  THE  PRESENT 

C PROGRAM  PERMITS  THE  USE  OF  TRACER  GAS  MEASUREMENTS  IN 

C DETERMINING  OFF-CYCLE  LOSSES.  AND  HANDLES  HEATING  EOUIP- 

C MENT  WITH  SINGLE.  TWO-STAGE.  OR  STEP-MODULATED  INPUT. 

C 

DIMENSION  TITLE(20.2) . CONC ( 20 ) . TS ( 20 ) . CONCM( 20 ) .TSM(20) 
DATA  PHI.REFTRM  /0.7.70./ 

C 

C INPUT  NUMBER  OF  DATA  SETS 
C 

WRITE  (6.500) 

READ  (5.*)  NRUN 
DO  310  1=1 .NRUN 


18 
1 9 
20 
21 
22 

23 

24 


C 

C INITIALIZE  COUNTER  AND  HIGH-FIRE  VARIABLES 
C 

NPASS=1 

QS0FFH=0. 

QICFFH=0. 

0S0FFM=0. 


25 


0I0FFM=0. 


26 

27 

28 

29 

30 

31 

32 


QINH=0. 

EFYSSH=0. 

EFFYUH=0. 

QI0NH=0. 

C 

C I/O  SYSTEM  PARAMETERS 

C 


33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 


READ  (5.510)  TITLE 

WRITE  (6.520)  ( ( T I T LE ( 1 1 . J J ) . 1 1 = 1 . 20 ) . d J=  1 . 2 ) 
READ  (5.+)  IFB. INST . ICTRL. LOSOPT 
IF  ( IFB. EO. 1 ) WRITE  (6.530) 

IF  (IFB.EQ.2)  WRITE  (6.540) 

IF  (1FB.E0.3)  WRITE  (6.550) 

IF  (INST.E0.1)  WRITE  (6.560) 

IF  (INST.E0.2)  WRITE  (6.570) 

IF  (ICTRL. EO.O)  WRITE  (6,580) 

IF  (ICTRL. E0.1)  WRITE  (6.590) 

IF  (ICTRL. EO-2)  WRITE  (6.600) 

IF  (LOSOPT. EO.O)  WRITE  (6.610) 

IF  ( LOSOPT . EO. 1 ) WRITE  (6,620) 

C 

C I/O  DATA  FOR  MAXIMUM  FIRING  RATE  AND  CALL  SENLCS 

C 


49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 


READ  (5.*)  NSYS,IFUEL.HHV.0IN.0P.P£,BE.XC02S 

WRITE  (6.630)  NSYS . I FUE L , HHV . QI N . QP . PE . BE . XC02S 

READ  (5.*)  TSSSX.XC02F.TFSS.TF0N1 .TF0N2.TF0FF3.TF0FF4. TF0FF5 

WRITE  (6,640)  TSSSX . XC02F . T FSS . T FON 1 . TF0N2 . TF0FF3 . TF0FF4 . TF0FF5 

READ  (5.*)  TRA.Od.SFR.DF.DS.Y 

CALL  SENLOS  ( I FUE L . NSYS . XC02S . TSSSX . XC02F . T FSS . HHVA . AFR . OL . RT , 

2 OSSS.EFFYSS.TRA. IFB.SFR) 

WRITE  (3.650)  TRA.Od.SFR.DF.DS.Y 

READ  (5.*)  TON. TOFF. TOA 

WRITE  (6.660)  TON. TOFF. TOA 

PF=QP/QIN 

Cd=3.3 

IF  (IFB.EQ.2)  Cd=4.7 
IF  (INST.E0.1)  Cd=0. 

C 

C I/O  data  FOR  MINIMUM  FIRING  RATE.  CALL  SENLOS  AND  WEIGHT 

C 

IF  (ICTRL. EO.O)  GO  TO  10 
NPASS=2 

READ  (5.*)  0INM.XC02SM 
WRITE  (6.670)  0INM.XC02SM 

READ  (5, * ) TSSSXM. XC02FM, TFSSM.TFONIM, TF0N2M. TF0F3M. TF0F4M, TF0F5M 
WRITE  ( 6 . 640 ) TSSSXM . XC02FM . TFSSM . T F0N1 M . TF0N2M . TF0F3M . T F0F4M . 

2 TFCF5M 

RE/D  (5.*)  TRAM.SFRM.DFM.DSM 

CALL  SENLOS  ( I FUE L . NS V S . XC02SM . TSSSXM . XC02FM , TFSSM . HHVA . AFR . QL . 

2 RTM.OSSSM. EFYS3M . TRAM. IFB.SFRM) 
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WRITE  (6,680)  TRAM . SFRM . DFM . DSM 

CALL  WEIGHT  ( 0 1 N . E FF YSS . 0 1 NM . EF YSSM . PLOW . TOA LO . TOAH I ) 

RHIGH=1 . -RLOW 

WRITE  (6,690)  R LOW , RH IGH , TOALO , TOAH I 
TOA=TOAHI 
C 

C I/O  TRACER  GAS  DATA  AND  CALCULATED  OFF-CYCLE  LOSS  IF  L0S0PT=1 

C 

10  IF  (LOSOPT.EO.O)  GO  TO  20 

READ  (5,*)  (CONC(d) ,TS( J) , J=1 ,20) ,CONCI ,VTTT,TROOM 

CALL  OFFLOS  ( CONC . TS . CONC I , VTTT , TROCM , TOA , 0 1 N . OSOF FH . 0 1 OF FH ) 

WRITE  (6,700)  (CONC( J ),d=1.10),(TS(d),d  = 1,10), (CONC( J) ,d=11 ,20) , 

2 (TS(J),d=11.20) .CONC I . VTTT ,TR00M,0S0FFH,0I0FFH 
IF  (ICTRL.EO.O)  GO  TO  20 

READ  (5,*)  (CONCM(d) ,TSM( J) ,d=1 ,20) ,CONCIM,VTTTM,TROOMM 
CALL  OFFLOS  ( CONCM , TSM , CONC IM , VTTTM , TROOMM , TOA LO , 01 NM . QSOFFM , 

2 OIOFFM) 

WRITE  (6,710)  (CONCM( J) ,d=1 ,10) , (TSM( J) ,J=1 . 10) . (CONCM(d) ,d=11 ,20) 
2 , (TSM( J) ,d=1 1 .20) . CONC IM. VTTTM, TROOMM. QSOFFM, OIOFFM 
20  WRITE  (6.720) 

C 

C LOOP  THRU  B0IL6  WORKING  CALCULATIONS  NPASS  TIMES 

C 

DO  270  J=1 .NPASS 
IF  ( J . EO. 1 ) GO  TO  30 
C 

C HOLD  NECESSARY  VARIABLES  FOR  MODULATING  AND  TWO  STAGE  CALCULATIONS 

C 

QINH=0IN 

EFYSSH=EFFYSS 

EFFYUH=EFFYU 

OIONH=OION 

C 

C REPLACE  MAXIMUM  FIRE  VARIABLES  WITH  MINIMUM  FIRE  VARIABLES 

C 

OIN=OINM 

TFSS=TFSSM 

TF0N1 =TF0N1M 

TF0N2=TF0N2M 

TF0FF3=TF0F3M 

TF0FF4=TF0F4M 

TF0FF5=TF0F5M 

TRA=TRAM 

SFR=SFRM 

DF=DFM 

DS=DSM 

RT=RTM 

OSSS=QSSSM 

EFFYSS=EFYSSM 

TOA=TOALO 

WRITE  (6,730) 

30  REFTOA=TOA 

WRITE  (6,740)  PF . HHVA . AFR . OL . CJ . RT . QSSS . EFFYSS 
C 

C ***  COLUMNS  31  THROUGH  43  IN  NBSIR  78-1543  ♦♦♦ 

C 

TSSS=(TFSS-TRA)/SFR+TRA 
IF  ( IFB. EQ.2)  GO  TO  40 
Cl  =2. 

C2=0. 5 
C3=7.5 
C4=1 .5 
GO  TO  50 
40  C1=4.5 

C2=1  . 

C3=18.75 
C4=3. 75 

50  IF  (TFSS.E0.TF0N1 . AND.TFSS.E0.TFON2)  GO  TO  60 

TA0N=C1/AL0G( ( TFSS-TF0N1 ) / ( T FSS - T F0N2 ) ) 

ZETF0X=( TFSS-TF0N1  ) * EXP ( C2/TA0N ) 

GO  TO  70 
60  TA0N=0. 

ZETF0X=0. 

70  CONTINUE 

TA0FF=C3/AL0G( ( T F 0FF3 - T FOF F5 ) / ( T FOF F4 - T F0FF5 ) ) 
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SIF0X=(TF0FF3-TF0FF5)*EXP(C4/TA0FF) 

SIF1X=TF0FF5-TRA 
IF  (NSYS.GT.8)  GO  TO  90 

IF  (NSYS.GT.4. AND.DS. LE. (DF/SFR) ) GO  TO  80 
DSF=DF/(SFR»DS) 

SISIX=DSF*SIFIX 
SISOX=DSF*SIFOX 
GO  TO  100 
80  SISIX=SIFIX 

SISOX=SI FOX 
GO  TO  100 
90  SISIX=0. 

SISOX=0. 

100  CS=1 .+(REFTRM-REFTOA)*EFFYSS/( 100. ‘(TFSS-REFTRM) ) 

IF  (NSYS.LT.9)  CS=0. 

CS0N=24.»(1 .+RT»AFR)/HHVA 
IF  (NSYS.GT.8)  GO  TO  120 
IF  (NSYS.GT.4)  GO  TO  110 

CSOFF=DF*CS0N*(TFSS+460.+REFTRM-TRA)**1 . 1 9/ ( T FSS - TRA ) * *0 . 56 
GO  TO  130 

110  CS0FF  = DS'»SFR»CS0N»(TSSS+460.+REFTRM-TRA)**1  . 1 9/ ( TSSS- T RA+REFTRM- 

2 REFT0A)**0.56 
GO  TO  130 

120  CSOFF  = DF*CSON* ( T FSS+460 . +REFTRM- TRA ) ♦ * 1 . 1 9/ ( T FSS - TRA+R EFTRM- 

2 REFTOA) »*0. 56 

130  IF  (NSYS.GT.8)  GO  TO  140 
CION=PHI*SFR»CSON 

CIOFF=DS*CION* ( TSSS-TRA+530. ) **1 . 1 9/ ( TSSS - TRA+RE FTRM- REFTOA ) • *0 . 56 
GO  TO  150 
140  CI0N=0. 

CI0FF=0. 

150  WRITE  (6.750)  TSSS . TAON . 2ETF0X . TAOFF . SI  FOX . S I FIX . SISIX . SI  SOX 
WRITE  (6.760)  CS . CSON . CSOF F . C I ON , C lOFF 
C 

C COLUMNS  44  THROUGH  53  IN  NBSIR  78-1543 

C 

IF  ( (HHV/HHVA) . LE.0.95)  WRITE  (6.770) 

IF  ( (HHV'/HHVA)  .GE.  1 .05)  WRITE  (6,780) 

IF  (TFSS.EO.TFON1 . AND.TFSS.E0.TF0N2)  GO  TO  160 
TTON=TON/TAON 
GO  TO  170 
160  TT0N=10.**20 

170  TT0FF=T0FF/TA0FF 

FON=SIFOX*EXP( -TTOFF)/( TFSS-TF0FF5) 

F0FF=2ETF0X*EXP( - TTON ) / ( TFSS- T F0FF5 ) 

F0N0F=1 . - FON'FCFF 
FF0N=(1 . -F0N)/F0N0F 
FFOFF=( 1 . -FOFF )/FONOF 

C CII=.90  IF  no  EQUIPPED  UNIT  ♦** 

IF  (OP.LT.0.1)  FFOFF=FFOFF* .90 

IF  (NSYS.GT.8)  GO  TO  180 

2ETF0=FF0N*ZETF0X 

SIFO=FFOFF*SIFOX 

SIFUSIFIX 

SISO=FFOFF*SISOX 

SISI=SISIX 

GO  TO  190 

180  2ETF0=CS*FF0N-2ETF0X 
SIF0=1 . 22»SI FOX^FFOFF 
SIFI=1 .22»SIFIX 
SIS0=0. 

SISI=0. 

190  WRITE  (6.790)  TOA . TON . TOF F . TTON , TTOFF . 2ETF0 . S I FO  . S I FI 
C 

C COLUMNS  54  THROUGH  59  IN  NBSIR  7B-1543 

C 

IF  (NSYS.GT.8)  GO  TO  210 
IF  (NSYS.GT.4)  GO  TO  200 
N = 4 

CALL  FUNT4  ( N . SI FO . TTOFF . REFTOA . REFTRM . F3 . F4 ) 

F5  = 0. 

F6  = 0. 

N = 8 

CALL  FUNT4  ( N . S I SO . TTOFF . REFTOA . REFT RM . F7 . F8 ) 

OSON=OSSS-CSON»2ETFO* ( 1 . -EXP( -TTON) )/TTON 
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0S0FF  = CS0FF* ( F3+S I F I * F4 ) *T0FF/T0N 
0ION=CION*( REFTRM-REFTOA) 

OIOFF=CIOFF*(REFTRM-REFTOA)*( F7+SISI*F8)*T0FF/T0N 
GO  TO  220 
200  CONTINUE 
F3  = 0. 

F4  = 0. 

N = 6 

CALL  FUNT4  ( N . S I SO . TTOFF . RE FTOA , REFTRM , F5 . F6 ) 

N = 8 

CALL  FUNT4  ( N . S I SO . TTOF F , REFTOA . REFTRM . F7 . F8 ) 
QSON=0SSS-CSON*2ETFO* ( 1 . -£XP( -TTON) )/TTON 
QSOFF=CSOFF* ( F5+S1 SI » F6 ) *TOFF/TON 
0ION=CI ON* {REFTRM- RE FTOA) 

OIOFF=CIOFF*(REFTRM-REFTOA)*( F7+S IS I* F8)* TOFF/TON 
GO  TO  220 
210  CONTINUE 
F3  = 0. 

F4  = 0. 

F7  = 0. 

F8  = 0. 

N = 6 

CALL  FUNT4  ( N . S I FO . JTOF F . REFTOA . REFTRM , F5 , F6 ) 

OSCN  = CS*QSSS-CSON*ZETFO*( 1 . - EXP ( -TTON) )/TTON 

QSOFF=CSOFF*( F5+SI FI*F6)*T0FF/T0N 

0ION=0. 

0IOFF=0. 

220  CONTINUE 
C 

C USE  MEASURED  VALUES  FOR  OFF-CYCLE  LOSS  IF  L0S0PT=1 

C 

IF  (LOSOPT.EO.O)  GO  TO  240 
IF  (J.E0.2)  GO  TO  230 
QSOFF=OSOFFH 
OIOFF=OIOFFH 
GO  TO  240 

230  QSOFF=OSOFFM 
OIOFF=OIOFFM 

240  IF  (NSYS.GT.8)  QI0FF=0. 

C 

C FINISH  UP  B0IL6  CALCULATIONS 

C 

IF  (INST.E0.2)  GO  TO  250 

EFFYU=100. -QL-TON* (OSON+OSOFF+OION+QIOFF)/( TON+PF*TOFF ) 

GO  TO  260 

250  EFFYU=100. -OL-CJ*Od- TON* (OSON+OSOFF )/(TON+PF* TOFF) 

260  WRITE  (6,800)  S I SO , S I S I . F3 . F4 . F5 . F6 . F7 . F8 
WRITE  (6.810)  OSON.OSOFF.OION.QIOFF.EFFYU 
IF  (ICTRL.GE. 1 ) GO  TO  270 

EFFYA=EFFYU*EFFYSS*4400/( EFFYSS*4400+2. 5*EFFYU*PF*4600) 

WRITE  (6.820)  EFFYA 
GO  TO  300 
270  CONTINUE 
C 

C CALCULATE  WEIGHTED  EFFICIENCIES  FOR  MODULATING  AND  TWO-STAGE 

C 

IF  (ICTRL.EQ.2)  GO  TO  280 
EFYSSA=( EFYSSH+EFYSSM)/2. 

QIONA=(OIONH+OION)/2 . 

EFYU=EFYSSA-OIONA 
EFYUW=EFYU*RHIGH+EFFYU*RLOW 
EFYSSW=EFYSSA*RHIGH+EFYSSM*RLOW 
GO  TO  290 

280  EFYUW=EFFYUH*RHIGH+EFFYU*RLOW 

EFYSSW=EFYSSH*RHIGH+EFYSSM*RLOW 

EFYU=EFFYUH 

290  EFFYAW=EFYUW*EFYSSW*4400/( EFYSSW*4400+2.5*EFYUW*PF*4600) 
WRITE  (6,830)  E FYU . EFFYU . EFYUW 
WRITE  (6.840)  EFYSSH . EFYSSM . EFYSSW 
WRITE  (6.850)  EFFYAW 
C 

C RELOAD  HIGH-FIRE  VARIABLES 

C 

0IN=01NH 

EFFYSS=EFYSSH 


Figure  C2  (cont) 


46 


EFFYU=EFYUW 


301 

302 

C 

303 

C 

304 

c 

305 

300 

306 

307 

308 

309 

310 

31  1 

312 

313 

314 

310 

315 

316 

31  7 

C 

318 

C 

319 

C 

320 

500 

321 

510 

322 

520 

323 

. 530 

324 

540 

325 

550 

326 

560 

327 

570 

328 

580 

329 

590 

330 

600 

331 

610 

332 

620 

333 

630 

334 

335 

336 

640 

337 

338 

339 

650 

340 

341 

660 

342 

670 

343 

344 

345 

346 

680 

347 

348 

690 

349 

350 

351 

352 

700 

353 

354 

355 

356 

357 

358 

359 

360 

361 

710 

362 

383 

364 

365 

366 

367 

368 

369 

370 

720 

371 

730 

372 

740 

373 

374 

750 

375 

CALCULATE  ANNUAL  FUEL  REQUIREMENTS 


A= 100000./ (34 1300. * ( PE+Y* BE )+ ( QIN -OP ) ♦ EFFYU ) 
B=2 . *A*QP*EFF YU/1  00000. 

IF  (OP. LT .0. 1 ) B=0. 

C=0.3B 

DHR=EFFYSS*QIN/ 100000. 

B0H=2080. *( A*C+DHR-B) 

AFUEL=(OIN-QP) *B0H+8760.*0P 
AELEC=(PE+Y*BE)’*BOH 
WRITE  (6.860)  AFUEL.AELEC 
CONTINUE 
WRITE  (6,870) 

STOP 


FORMAT  STATEMENTS 


FORMAT  (/5X.31HINPUT  DATA  --  ©ADD  F I LE . ELEMENT/ ) 

FORMAT  (20A4) 

FORMAT  (//2(5X,20A4/) ) 

FORMAT  (/5X.7HFURNACE) 

FORMAT  (/5X.6HB0ILER) 

FORMAT  (/5X. 13HVENTED  HEATER) 

FORMAT  (5X. 16HINSTALLED  INDOOR) 

FORMAT  ( 5X.  1 7HINSTALLED  OUTDOOR) 

FORMAT  (5X. IBHSINGLE  STAGE  INPUT) 

FORMAT  (5X.20HSTEP  MODULATED  INPUT) 

FORMAT  (5X.15HTW0  STAGE  INPUT) 

FORMAT  ( 5X . 30HASS1GNED  VALUES  FOR  OFF-LOSSES) 

FORMAT  (5X.30HMEASURED  VALUES  FOR  OFF-LOSSES) 

FORMAT  (///5X. 13H1NPUT  VALUES  //5X . 6H1 ) NS YS4X . 7H2 ) I FUE L3X . 

2 5H3)HHV5X.5H4)0IN5X.4H5)0P6X,4H6)PE6X.4H7)BE6X.7H8)XC02S/I7. 110, 

3 6X.8(1PE10.2)  ) 

FORMAT  (/5X .6H9)TSSS4X .8H10)XC02F2X .7H1 1 ) T FSS3X . 8H 1 2 ) T FON 1 2X . 

2 8H13)TF0N22X. 9H14)TF0FF31 X.9H15)TF0FF41X.9H1 6)TF0FF5/3X. 

3 8( 1PE10.2)  ) 

FORMAT  ( /5X . 6H 1 7 ) T R A4X . 5H 1 8 ) 0U5X . 6H 1 9 ) S/F4X , 5H20 ) D F5X . 5H2 1 ) DS . 5X . 

2 8H22)Y  /3X,8( 1 PE10.2) ) 

FORMAT  (/5X.4HT0N=F6.2.4X.5HT0FF=F6.2,5X,4HT0A=F5. 1 ) 

FORMAT  (///5X.47HADDITI0NAL  INPUT  VALUES  FOR  MINIMUM  FIRING  RATE// 

2 5X.6H1 )NSYS4X . 7H2 ) I FUE L3X . 5H3 ) HHV5X , 5H4 )QI N5X . 4H5 ) QP6X , 4H5 ) PE6X . 

3 4H7)BE6X.7H8)XC02S/3(7X,3H--- ) . 3X . 1 PEI  0 . 2 . 4X . 3H - - - . 2 ( 7X . 3H - - - ) . 

4 3X.1PE10.2) 

FORMAT  ( /5X . 6H1 7 ) TRA4X . 5H1 8 )0J5X . 6H1 9 ) S/F4X . 5H20 ) DF5X . 5H2 1 )DS . 5X , 

2 8H22)Y  /3X. 1PE10.2.4X.3H--- ,3X,3( 1PE10.2) .4X.3H--- ) 


FORMAT  (//5X.49HFRACTI0N  OF  TIME  IN  LOW  FIRE/CYCLING  MODE  = . 

2 F5.3/5X.49HFRACTI0N  OF  TIME  IN  HIGH  F I RE/NON - C YC L I NG  MODE  = F5.3/ 

3 5X.40H0UTD00R  AIR  TEMP  FOR  CYCLING  MODE  = F6.2/5X, 

4 40H0UTD00R  AIR  TEMP  FOR  NON-CYCLING  MODE  = F6.2) 

FORMAT  (///5X.22HFLUE  LOSS  MEASUREMENTS// 


63H  TIME 
13H  8-9 

63H  TIME 
13H  18-19 
26HINITIAL 


0-1  1-2 
9- 10/6H  CONC 
10-11  11-12 
19-20/6H  CONC 


2-3  3-4  4-5 

10( IX. F6. 1 )/6H  TEMP 
12-13  13-14  14-15 

10( IX. F6. 1 )/6H  TEMP 


5-6  6-7  7-8 

10( 1X.F6. 1 )// 

15-16  16-17  17-18 

10( IX. F6. 1 )//5X. 


CONC( ENTRATION)  = 


1PE10.2/5X. 

27HMETER  VOLUME/BUBBLE  TIME  = 1PE10.2/5X. 

24HTEST  ROOM  TEMPERATURE  = 1PE10.2//5X. 

17HQS0FF  MEASURED  = 1 PE  1 0 . 2/5X . 1 7HQI0FF  MEASURED  = 
FORMAT  (//5X, 19HMINIMUM  FIRE  VALUES// 


63H  TIME 
13H  8-9 

63H  TIME 
13H  18-19 
26HINITIAL 


0-1  1-2  2-3 

9- 10/6H  CONC  10( IX, 
10-11  11-12  12-13 

19-20/6H  CONC  10( IX, 
C0NC( ENTRATION) 


3-4  4-5 

F6. 1 )/6H  TEMP 
13-14  14-15 

F6, 1 )/6H  TEMP 


1 PE10.2) 
6-7 

F6. 1 )// 
16-17 


7-8 


5-6 
10(  IX 
15-16  16-17  17-18 

10(1X.F6.1 )//5X, 


1 PE10. 2/5X, 

27HMETER  VOLUME/BUBELE  TIME  = 1PE10.2/5X. 

24HTEST  ROOM  TEMPERATURE  = 1PE10.2//5X. 

17H0S0FF  MEASURED  = 1 PEI  0 . 2/5X , 1 7HQI0FF  MEASURED  = 1PE10.2) 

FORMAT  ( ///5X , 1 8HCALCULATED  VALUES-/) 

FORMAT  (///5X .40HCALCULATED  VALUES  --  MINIMUM  FIRING  RATE/) 

FORMAT  (/5X.5H23)PF5X,7H24 )HHVA3X ,6H25)A/F4X.5H26)0L5X,5H27)CJ5X. 

2 5H23)RT5X.7H29)0SSS3X.9H30)EFFYSS/3X.8(1PE10.2) ) 

FORMAT  ( /5X . 7H3 1 ) TSSS3X . 8H32 ) TAU0N2X . 9H33 ) ZET  FOX  1 X . 9H34 ) T AUOFF , 1 X . 
2 9H35)PSIF0X1X ,9H36)PSIFIX1X.9H37)PSISIX1X,9H38)PSIS0X/3X . 
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IS*B0IL6A( 

1 ) .1^ 

1 

2 

C 

3 

C 

4 

C 

5 

C 

6 

c 

7 

c 

8 

9 

10 

1 1 

1 2 

c 
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1 9 
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33 

7H47)TT0N3X. 

IPSIFI/3X, 


3 1PE10.2) ) 

FOR'.UT  ( /5X  . 5H39)CS5X  . 7H40  ) CSON3X  . 8H4 1 )CSOFF2X  . 7H42  )CI0N3X  , 

2 8H43)CI0FF/3X.5( 1 PE10.2) ) 

FORMAT  (///5X.46H***  WARN  I NG- HEAT  I NG  VALUE  OF  TEST  FUEL  IS  TOO 
2 7HLOW  ***///) 

FORMAT  (///5X.46H***  WARN  I NG- HEAT  I NG  VALUE  OF  TEST  FUEL  IS  TOO 
2 8HHIGH  ***///) 

FORMAT  ( /5X . 6H44 ) T0A4X . 6H45 ) TON4X . 7H46  ) TOFF3X , 

2 8H48)TT0FF.2X,8H49)ZETF02X.8H50)PS1F02X,8H51 : 

3 8( 1PE10.2)  ) 

FORMAT  ( /5X . 8H52 ) PSIS02X . 8H53)PSISI2X , 5H54 ) F35X . 5H55) F45X . 

2 5H56)F55X.5H57)F65X.5H58)F75X,5H59)F8/3X.8( 1 PE  10. 21) 

FORMAT  ( /5X . 7H50 ) QS0N3X . 8H6 1 )0S0FF2X . 7H62)0I0N3X . 8H63)QI0FF2X . 

2 8H64)EFFYU/3X,5( 1 PE10.2) ) 

FORMAT  (///5X.25HANNUAL  FUEL  EFFICIENCY  = F6.2) 

FORrMT  ( ///5X . 35HPART - LOAD  EFFICIENCY  --  HIGH  FIRE  = F6.2/5X, 

2 35HPART-L0AD  EFFICIENCY  --  LOW  FIRE  = F6.2/5X, 

3 40HWEIGHTED  AVERAGE  PART-LOAD  EFFICIENCY  = F6.2) 

FORMAT  {/5X.39HSTEADY-STATE  EFFICIENCY  --  HIGH  FIRE  = F6.2/5X, 

2 33HSTEADY-STATE  EFFICIENCY  --  LOW  FIRE  = F6.2/5X. 

3 43HWEIGHTED  AVERAGE  STEADY-STATE  EFFICIENCY  = F6.2) 

FORMAT  (/5X.42HWEIGHTED  AVERAGE  ANNUAL  FUEL  EFFICIENCY  = F6.2) 
FORMAT  (5X.33HANNUAL  FUEL  ENERGY  CONSUMPTION  = 1 PEI  0 . 3 . 1 X . 3HBTU/ 
2 5X.39HANNUAL  ELECTRICAL  ENERGY  CONSUMPTION  = 1 PEI  0 . 3 . 1 X . 3HBTU ) 
FORMAT  (///5X. 16HPR0GRAM  COMPLETE//) 

END 


Subroutine  I-7EIGHT 


SUBROUTINE  WE  I GHT ( 01 N . EFFYSS . 01 NM . E FYSSM , R LOW . TCALO . TO AH  I ) 

SUBROUTINE  TO  DETERMINE  THE  FRACTION  OF  HEATING  SEASON  IN 
MINIMUM  FIRE/CYCLING  MODE  AND  IN  MAXIMUM  FiRE/NON-CYCLING 
MODE, AS  WELL  AS  THE  CORRESPONDING  OUTDOOR  AIR  TEMPERATURE 
FOR  EACH  MODE. I .E. .TOALO  AND  TOAHI 

DIMENSION  RL0(9)  . TOALO)  ,T0AH(9) 

DATA  RLO/0.01 .0.09.0.21.0.33.0.45.0.57.0.69.0.81 .0.89/ 

DATA  TOAL/59.4.57.9.55.2.53.0.51 .1.49.5.47,7.45.9.45.0/ 
DATA  TOAH/42. 1,40.6.38.2.36.1.33.9.31.6.28.6.24.3.21.7/ 


HIMAX=OIN* EFFYSS/ 100. 

himin*oinm*efyssm/ioo. 

RATIC=HIMIN/HIMAX 
N = 0 

ASSIGN  VALUES  FOR  RLOW.  TOALO.  AND  TOAHI 


IF(RATIO. LE.O. 
IF(RATI0. LE.O, 
IF(RATI0. LE.O, 
IF(RATIO. LE.O. 
IF(RATIO. LE.O. 
1F(RATI0. LE.O, 
IF(RATIO. LE.O. 
IF(RATI0. LE.O. 
IF(RATIO.GT.O. 
RLOW=RLO(N) 
TOALO=TOAL(N) 
TOAHI=TOAH(N) 
RETURN 
END 


15)  N=1 

26. AND. RATIO. GT.O. 15)  N=2 
35. AND. RATIO. GT. 0.26)  N=3 
43. AND. RATIO. GT.O. 35)  N=4 
50. AND. RATIO. GT.O. 43)  N=5 
57. AND. RATIO. GT .0. 50)  N=6 
66. AND. RATIO. GT.O. 57)  N=7 
79. AND. RATIO. GT.O. 66)  N=8 
79)  N=9 


Figure  C2  (cont) 
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Subroutine  SENLOS 


HTRS*BOI L6A( 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 
1 7 
18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


1 ) .SENLOS( 1 ) 

SUBROUTINE  SENLOS ( I FUE L . NSYS . XC02S . TSSSX . XC02 F . T FSS , HHVA , APR  . OL  , RT 
* .QSSS.EFFYSS.TRA.  1 FB.SFR) 

C ***CALCULATI0N  OF  HHVA  AFR  OL  RT  QSSS  EFFYSS  *+* 

C ***APRIL  1978+**  FOR  USE  WITH  NBSFBS  5 ♦** 

DIMENSION  HHV(6) .AF(6)  ,0(6) ,ART(6) ,BRT(6) . CA ( 5 ) . CF ( 6 . 5 ) 

DATA  (HHV(J) ,d=1 ,6)/19800. . 19500. .20120. . 18500. .21500. .20890./ 

DATA  (AF(J).J=1,6)/14.56.14.49.14.45.11.81.15.58, 15.36/ 

DATA  (0( J) .0=1 .6)/6.55.6.50.9.55. 10. 14,7.99.7.79/ 

DATA  (ART(J) , J=1 .6)/. 0679. .06668. .09194, .09646. .0841 . .0808/ 

DATA  (GRT(d) ,d=1.6)/14.22.14.34.10.96,10.10, 12.60.12.93/ 

DATA  (CA(d) .d=1 .5)/2.5462121E-1 . - 3 . 02601 26E - 5 . 2 . 760857 1 E - 8 , -7.4253 
♦321E-12.6.4307377E-16/ 

DATA  (CF( 1 ,K) ,K  = 1 . 5)/2 . 441 6834E-01 . 3 . 37 1 1 449E - 6 . 8 . 8906305E - 9 . -1.36 
* 1 901 SE- 12,-1  . 43674 10E- 16/ 

DATA  (CF(2.K) ,K=1  . 5 ) /2 . 4361  1 63E- 1 . 3 . 6702686E - 6 , 8 . 7098897E - 9 . - 1 .309 
*437SE-12. -1 .5029209E- 16/ 

DATA  (CF(3.K)  , K=1  , 5 ) /2 . 5949478E- 01 . - 4 . 9475802E - 06 . 1 .3885838E-8. -2. 
*C059994E- 1 2.3. 7682444E - 1 7/ 

DATA  (CF(4.K) . K=1  . 5 ) /2 . 6598442E- 1 . - 7 . 7561 435E -6 . 1 . 5833852E - 8 . - 3 . 41 
♦94210E-12.1 .2158977E-16/ 

DATA  (CF(5,K)  ,K=1  , 5 ) /2 . 51 63639E- 1 . - 6 . 4 1 44604E - 7 . 1 . 1 31 5073E-8 . -2 . 06 
♦56792E-12.-5.4897330E-17/ 

DATA  (CF(6.K) ,K=1  . 5 ) /2 . 501  1 247E - 1 ,1 . 7737005E - 7 . 1 . 0820337E-8 . - 1 .922 
*0641 E- 12. -7.3013274E- 17/ 

HHVA=KHV( IFUEL) 

AFR=AF( IFUEL) 

QL=0( IFUEL) 

RT=(BRT( IFUEL)/XC02F)+ART{ IFUEL) 

XC02=XC02F 

TSS=TFSS+459.69 

C ♦**  SFR  TEST  ONLY  FOR  VENTED  HEATERS  WITH  DRAFT  HOODS  *** 

IF(IFB.NE.3.0R.XC02S. LT.O. 1 ) GO  TO  30 
RTS=  (BRT( IFUEL)/XC02S)+  ART(IFUEL) 

CALSFR  = 1 .3*(RTS/RT) 

SFR  = CALSFR 

30  IF(XC02S. LT .0. 1 .OR. TSSSX. LT.O. 1 ) GO  TO  50 

XC02=XC02S 
TSS=TSSSX+459.69 
50  CONTINUE 

RTX  = ART(  IFUEL)  + (BRT(  IFUED/XC02) 

Qr  = 0. 
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0A  = 0. 


43 

44 

45 

46 

47 

48 

49 

50 


DO  100  1=1.5 

QF  = QF+CF( IFUEL. I ) * ( TSS* * I - ( TRA+459 . 69 ) * * I ) 
100  OA=QA+CA( I)*(TSS**I-( TRA+459 . 69 )♦* I ) 
QSSS=(1 .+AFR)*0F+(RTX-1 . )*AFR*QA 
0SSS=100. *0SSS/HHVA 
EFFYSS=100. -QL-OSSS 
RETURN 
END 


Subroutine  OFFLOS 


HTRS+BOI L6A( 1 ) .0FFL0S(2) 

1 SUBROUTINE  OFF  LOS ( CONC . TS . CONC I . VTTT , TROOM , TOA . QIN . QSOFF , OIOFF ) 

2 C 

3 C SUBROUTINE  TO  CALCULATE  OFF-CYCLE  FLUE  LOSSES  BASED  ON  MEASURED 

4 C VALUES  OF  TRACER  GAS  CONCENTRATION  AND  STACK  TEMPERATURE 

5 C 


7 

8 
9 

10 
1 1 
1 2 
1 3 
1 4 
1 5 
16 
1 7 
18 


DIMENSION  C0NC(20) ,TS(20) 

0FFLS=0.0 
0FFLI=0.0 
DO  10  J=1 ,20 

R,ATI0=(C0NCI  -CONC(U)  )/CONC(d) 
r LCL'^O.  01  * VTTT’ RAT  10/60. 
0FFLS=0FFLS+0 . 24* F LOW* ( TS ( J )- TROOM ) 
0FFLI=0FFLI+0. 24*FL0W*( 70.0-T0A) 

10  CONTINUE 

QS0FF=3C0. *0FFLS/QIN 
0IOFF=30O,*1 .3*0. 7*OFFLI/0IN 
RETURN 
END 


Figure  C2  (cont) 
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Subroutine  FUNT4 


HTRS*B0IL6A( 1 ) . FUNT4( 1 ) 


1 

2 

3 

4 
B 
6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 
1 9 
20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


100 


110 

120 


SUBROUTINE  FUNT4(N.SI . TTOFF . REFTOA . REFTRM, FI , F J ) 

T1=REFTRM-REFT0A 

C1  = 1 . 

T2=REFTRM+460. 

IF(N.E0.4)  T1=0. 

IF(N.EQ.8)  C1=0.0 
X = 0. 

FI  =0. 

FU  =0. 

DX  =TTOFF  /500. 

FF1  1 = (SI  + T1 )**0.56*SI**C1/(SI+T2)**1 .19 

FF21  = (SI+T1+100.  )**0.56’»(SI  + 100.  ) ♦ *C1  / ( SI+T2+1 00 . )**1  . 19 

DO  100  1=1 .250 

X=X+DX 

XX=SI  *EXP(-X) 

FF12=(XX+T1 )**0.56*XX**C1/(XX+T2)**1 .19 

FF22=(XX+T1+100. ) * *0 . 56* ( XX+1 00 . ) • *C1 / ( XX+T2+ 1 00 . )**1 . 19 
X=X+DX 

XX=SI  *EXP(-X) 

FF13=(XX+T1 )**0.56*(XX)**C1/(XX+T2)**1 . 19 
FF23=(XX+T1+100. ) * *0 . 56* ( XX+1 00 . )**C1/(XX+T2+100. )**1 . 19 
FI  =FI  +(FF11+4.*FF12+FF13) 

Fd  =FJ  +(FF21+4.*FF22+FF23) 

FF1 1=FF13 
FF21=FF23 

IF( (ABS( Fd-FI  )). LE. 0.0000001 ) GO  TO  110 
Fd  =(Fd  -FI  )/( 100. *TTOFF  ) 

Fd  =DX  *Fd  /3. 

GO  TO  120 
Fd  =0. 

FI  =FI  *DX  /(3.*TT0FF  ) 

RETURN 

END 


Figure  C2  (cent) 
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Figure  C3 . Input  data  format  for  FBVH 


Line  1:  NRUN  (integer) 

NRUN  = number  of  sets  of  test  data  to  be  analyzed.  Lines  2 through  53 
should  be  repeated  NRVN  times  in  the  stored  data  file/element 

Line  2:  TITLE 

TITLE  = one  to  80  alphanumeric  characters  describing  the  test  conditions/data 


Line  3:  SUBTITLE 

SUBTITLE  = one  to  80  alphanumeric  characters  describing  the  test  conditions/data 


Line  4:  IFB,  INST,  ICTRL,  LOSOPT  (all  integers) 

IFB  = 1 for  furnaces 
= 2 for  boilers 
= 3 for  vented  heaters 

INST  = 1 for  vented  heating  equipment  installed  indoors 
= 2 for  vented  heating  equipment  installed  outdoors 
and  floor  furnaces 

ICTRL  = 0 for  single  stage  controls 

= 1 for  step  modulating  controls 
= 2 for  two  stage  controls 

LOSOPT  = 0 for  analytic  determination  of  off-cycle  stack  losses 
(method  used  in  References  1 and  4) 

= 1 for  determination  of  off-cycle  stack  losses  from  tracer 
gas  technique  (see  Reference  6) 

Line  5:  NSYS,  IFUEL,  HHV,  QP,  PE,  BE,  XC025 

Line  6:  TSSS,  XC02F,  TFSS,  TFONl,  TFQN2,  TF0FF3,  TF0FF4,  TF0FF5 
Line  7:  TRA,  Q J , SFR,  DF,  DS , Y 

Measured  quantities/assigned  values  defined  in  Reference  4,  evaluated 
at  maximum  firing  rate 

Line  8:  TON,  TOFF,  TOA 

TON  = 20  min 

TOFF  = 20  min  For  Vented  Heaters 
TOA  = 45°F  J 

See  Reference  4 for  values  for  other  types  of  heating  equipment. 
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OMIT  LINES  9-11  for  single  stage  controls  (ICTRL  = 0) 

Line  9:  QINM,  XC02SM 

Line  10:  TSSSM,  XC02FM,  TFSSM,  TFONIM,  TFON2M,  TFOF3M,  TFOF4M,  TFOF5M 

Line  11:  TRAM,  SFRM,  DFM,  DSM 

Measured  quantities /assigned  values  evaluated  at  minimum  firing  rate 
OMIT  LINES  12-53  for  Analytic  Determination  of  Off-Cycle  Losses  (LOSOPT  = 0) 

Lines  12-31-  CONC,  TS 

CONC  = measured  concentration  of  tracer  gas  in  sample  taken  from  stack,  PPM 
TS  = measured  stack  gas  temperature,  “F 

Twenty  sets  of  concentration  and  temperature  measurements  for  a cooldown  following 
steady  state  operation  at  maximum  fire.  For  more  details  on  tracer  measurements 
see  Reference  3. 

Line  32:  CONCI,  VTTT,  TROOM 

CONCI  = concentration  of  active  traver  gas  in  the  tracer  gas  in  PPM,  eg, , for 
a carbon  monoxide/nitrogen  tracer  gas  mixture  consisting  of  50.8 
percent  carbon  monoxide,  CT  = 508,000  PPM 

3 

VTTT  = volumetric  flow  rate  of  tracer  gas,  cm  /sec 

TROOM  = test  room  temperature,  **F 
Values  for  test  at  maximum  firing  rate. 

OMIT  LINES  33-53  For  Single  Stage  Controls  (ICTRL  = 0) 

Lines  33-52:  CONCM  TSM 

Twenty  sets  of  concentration  and  temperature  measurement  for  a cooldown  following 
steady  state  operation  at  minimum  fire. 

Line  53:  CONCIM,  VTTTM,  TROOMM 

Values  for  test  at  minimum  firing  rate. 


Figure  C3  (cont) 
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MAXIMUM  INPUT  MINIMUM  INPUT  MAXIMUM  INPUT 


Figure  C4.  Input  data  form  for  FBVH 


LINE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


f ~~  FBVH  DATA~  ~ | ^N  IK  | DATE 

NRUN 


TITLE 


SUBTITLE 


IFB 

INST 

ICTRL 

LOSOPT 

NSYS 

IFUEL 

HHV 

QIN 

QP 

PE 

BE 

XC02S 

TSSS 

XC02F 

TFSS 

TFONl 

TF0N2 

TF0FF3 

TF0FF4 

TF0FF5 

TRA 

QJ 

SFR 

DF 

DS 

Y 

TON 

TOFF 

TOA 

QINM 

XC02SM 

TSSSM 

XC02FM 

TFSSM 

TFONIM 

TF0N2M 

TF0F3M 

TF0F4M 

TF0F5M 

TRAM 

SFRM 

DFM 

DSM 

'12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


TIME 

CONC 

TS 

0-1 

1-2 

2-3 

3-4 

4-5 

5-6 

6-7 

7-8 

8-9 

9-10 

10-11 

11-12 

12-13 

13-14 

14-15 

15-16 

16-17 

17-18 

18-19 

19-20 

CONCI 

VTTT 

TROOM 

LINE 

’33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


TIME 

CONCM 

TSM 

0-1 

1-2 

2-3 

3-4 

4-5 

5-6 

6-7 

7-8 

8-9 

9-10 

10-11 

11-12 

12-13 

13-14 

14-15 

15-16 

16-17 

17-18 

18-19 

19-20 

CONCIM 

VTTTM 

TROOMM  ] 
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Figure  C5.  Program  output  for  a heater  with  step-modulating  control 
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FRACTION  OF  TIME  IN  LOW  FIRE/CYCLING  MODE  * .690 

FRACTION  OF  TIME  IN  HIGH  FIRE/NON-CYCLING  MODE  ■ .310 

OUTDOOR  AIR  TEMP  FOR  CYCLING  MODE  = 47.70 

OUTDOOR  AIR  TEMP  FOR  NON-CYCLING  MODE  = 28.60 
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Figure  C6.  Program  output  for  a heater  with  two-stage  control 
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fraction  of  time  in  low  fire/cycling  mode  = ,690 

FRACTION  OF  TIME  IN  HIGH  FIRE/NON-CYCLING  MODE  * .310 

OUTDOOR  AIR  TEMP  FOR  CYCLING  MODE  = 47.70 

OUTDOOR  AIR  TEMP  FOR  NON-CYCLING  MODE  = 23.60 
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